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Resume 

Notre objectif est de mettre en evidence l'instabilite spectrale de cer- 
tains operateurs differentiels non-autoadjoints, via l'etude de la croissance 
des normes des projecteurs spectraux II n associes a leurs valeurs propres. 
Nous nous interessons a certains oscillateurs anharmoniques 
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definis sur Z/ 2 (R). Nous etendons les resultats de [4] et [5] en donnant un 
developpement asymptotique de la norme des projecteurs spectraux as- 
socies aux grandes valeurs propres pour les operateurs -4(1, 9) et A(2k, 9), 
k>l. 

The purpose of this Note is to highlight the spectral instability of some 
non-selfadjoint differential operators, by studying the growth rate of the 
norms of the spectral projections H n associated with their eigenvalues. 
More precisely, we are concerned with some anharmonic oscillators 

with 
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defined on L (R). We get asymptotic expansions for the norm of the 
spectral projections associated with the large eigenvalues of .4(1, 9) and 
A(2k,9), k > 1, extending the results of 2] an< i El- 



1 Spectral instability and pseudospectra 

It is well-known that the spectral theorem implies some control of stability 
for the spectrum of selfadjoint operators : if A is a selfadjoint operator acting on 
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the Hilbcrt space H, the spectrum of its perturbations A + eB, with e > and 
any B € \\B\\ < 1, lies entirely inside an e-neighborhood of the spectrum 

a(A). In other words, the norm of the resolvent of A near the spectrum blows 
up like the inverse distance to the spectrum. It has also been known for several 
years (see for instance [13]) that such a behavior could not be expected in general 
in the case of non-selfadjoint operators. One can understand it thanks to the 
notion of e-pseudospectra of an operator A, defined as the family of sets o~ e (A), 
indexed by e > 0, 



<Ts{A) = U&p{A):\\{A-0- 1 \\> Hua 



(A). 

The link between spectral instability and pseudospectra appears more clearly in 
the following equivalent formulation, which is a weak version of the Roch and 
Silbermann theorem [11] : 



a e (A) = \J a(A + sB) 



e £(•«), 
II J» II < i 



(See also [H] and references therein). 

In the following, we deal with the instability indices associated with an isolated 
eigenvalue A £ a (A). The instability index associated with A is defined as 

«(A) = ||n(A)||, 

where 11(A) denotes the spectral projection associated with A. Of course k(A) > 
1 in any case, and re(A) = 1 when A is selfadjoint. These numbers k(A) are 
closely related to the size of e-pseudospectra around A. Indeed, if a* denotes 
the connected component of o~ £ (A) containing A, and if we assume for simplicity 
that <7g Da(A) = {A} and a* is bounded, then the perimeter \d<r^\ of satisfies 
(see i) 

\do$\ > 2tt£k(A). (1.1) 

In the finite dimensional setting at least, instability indices give a better descrip- 
tion of pseudospectra : if A € A4 n (C) is a diagonalisable matrix with distinct 
eigenvalues Ai, . . . , A„, Embree and Trefethen show fTS] that there exists £o > 
such that, for all e < Eq, 

|J D(X k ,e K (X k )+0(e 2 ))c<i £ (A)c (J D(X k , £K (X k )+0(e 2 )). (1.2) 

Xh&<r(A) A fc e<r(yl) 

In the case of an infinite dimensional space, the validity of this statement should 
be investigated, as well as the dependance on X k of the 0(e 2 ) terms. 
In the following, we study the instability indices of simple non-selfadjoint dif- 
ferential operators introduced by Davies in [4], for which the instability phe- 
nomenon described above will appear clearly. Let us define the anharmonic 
oscillators ^ 

A(m,6) = -— + e ie \x\ m (1.3) 

with 

... f (m + 2)7r (m + 2)7rl 

|*|<nnn — .-^ , (1.4) 
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defined on L 2 (K) in [¥| by taking the closure of the associated quadratic form 
defined on Cq°(K), which is sectorial if 9 satisfies (|1.4[) . According to 0], its 
spectrum consists of a sequence of discrete simple eigenvalues, denoted in nonde- 
creasing modulus order by A n = \ n (m,9), |A„| — > +00. The associated spectral 
projections are of rank 1, and E.-B. Davies showed in [4] that for all m G]0, +oo[ 
and 9 7^ satisfying (jl.4[) . for all a > 0, there exists N = N(m, 9, a) > such 
that the instability indices K n (m, 9) of A(m, 9) satisfy K n (m, 9) > n a for n> N. 
This statement has been improved in the case m — 2 of the harmonic oscillator 
(sometimes refered as the Davies operator), since E.-B. Davics and A. Kuijlaars 
showed [5] that K n (2,6) grows exponentially fast as n — > +00, with an explicitc 
rate c(9) : there exists an explicit c(9) such that 

Urn -logK„(2,0) = c(0). (1.5) 

n— >+oo n 

The purpose of this Note is to prove that this statement actually holds for the 
so-called complex Airy operator A(l, 9) and for the even anharmonic oscillators 
A(2k,6), k > 1. 
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2 Non-selfadjoint anharmonic oscillators 

We first deal with the complex Airy operator A(l,9) defined in (|1.3[) . We 
show that the corresponding instability indices K n (l,9) grow like in (|1.5[) as 
n — > +00. More precisely, we get asymptotic expansions in powers of as 
n — > +00. 



Theorem 2.1. LetO < \9\ < Sir/A. There exists a real sequence (&j(9))j>i such 
that the instability indices K n (l,9) of A(l,9) satisfy, as n — > +00, 

exp {-C{6){n- 1/2)) Kn (1, 6) = ^=2 f 1 + ^ «j(«)n _i j +0(n-°°), (2.1) 
where 

C{9) = nrnj 2 \ sin 9\ and K{9) 



1/4 ' 



2^3|sin6»| 7 
with 

I ' sin 2 (20/3) o cos(0/3)sin(2g/3) 

m e = H ttt; 2 — > 0. 

y sur 9 sin 9 

Sketch of the proof : Let us first recall that all the eigenvalues of A(m, 9), 
m G N, have associated spectral projections of rank 1, see Lemma 5 in 0]. Hence, 
one can easily check that [3] 

K n (ra,9) = ^^, (2.2) 
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where u n denotes an eigenfunction associated with the n-th eigenvalue of A(m, 9). 
We get rid of the singularity of the potential at x = by decomposing ,4(1,0) 
into its Dirichlet and Neumann realizations A D (1,8) and A (1,0) on M. + . Wc 
then compute their instability indices 

/ R+ ^(/i^+e^x) 2 ^ 

given by formula (|2.2p . where x H> Ai([in^ N + e l8 / 3 x) is the n-th eigenfunction of 
^4.- D / JV '(l 5 0), (rcsp. jt*^) being the n-th (negative) zero (resp. critical point) 
of the Airy function Ai (see [2], [8]). We estimate the numerator in 1(2.30 by 
using the well-known asymptotic expansion of at infinity in the complex 
plane (see [I]), and the Laplace method brings an exp(ce|/Xn ^J 3 / 2 ) term in 
Kn^ N (l,9), eg > 0. The integral in the denominator of (|2.3|) . after deformation 
of the path of integration by homotopy, is equal to 

Ai 2 (x)dx = Ai/ 2 (^/ N ) (2.4) 

(it is indeed straightforward, using Airy equation, to check that x i-> xAi 2 (x) — 
Ai' 2 (x) is a primitive for Ai 2 ). Hence the expansion of Ai'(-z) as z — > +oo, 
given in [I], provides an asymptotic expansion for (|2.4p in powers of \fin / N \~ 3 / 2 . 
The statement follows from the behavior of fin^ N as n — > +oo, since we have 
asymptotic expansions for (n — l/4) _2 / 3 |/z^| (resp. (n — 3/4) _2 / 3 |u^|) in powers 
of (n - l/4)~ 2 (resp. (n - 3/4)- 2 ). □ 

Notice that the exponential instability appears as soon as 9 ^ 0. 
Wc have a similar statement for even anharmonic oscillators : 

Theorem 2.2. Let k e N* and 9 be such that < \9\ < (fc +^ 1)7r . If n n (2k,9) 
denotes the n-th instability index of A(2k,9) = —-jbi + e ie x 2k , then there exist 
K(2k,0) > and a real sequence (C J (2k, 0))j>i such that 

e- Ck{e)n n n (2k,9) = K( - 2l !l g) ( 1 + Yc ] (2k,9)n-J I + ( 2 .5) 



as n —> +oo ; with 



etW = 2( t+1 )^r(y )yt ( 6 ) >0i (26) 

1 V2/J 



where 



tan(|fl|/(fc + l)) 

Vsin(fc|0|/(fc + l))+cos(fc|0|/(fc + l))tan(|0|/(fc + l))/ *• ' ; 



P*(0 = Im / (l-i 2fe ) 1/2 * 
Jo 



(2. 
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Sketch of the proof : We first perform an analytic dilation and a scale 
change to recover the semiclassical selfadjoint anharmonic oscillator Vh{2k) = 



i 2 -^ + x 2k - 1, with 



h = h n = \\ n {2k,9)\-^. 



The n-th instability index of A(2k,8) then writes 

Kn (2k, 9) = jRim f K n 2.9 

where iph solves Vh(2k)tph = 0, t/jh S L 2 (R) (sec (|2.2|l . after deformation of the 
integration path in the denominator). The complex WKB method (see [TU], |14) . 
[B]) and the analysis of the Stokes lines of Vh{2k) provide an asymptotic expan- 
sion of iphie 1 2{k+1) x) as h — > 0, which enables us to determine the asymptotic 
behaviour of the numerator in (|2.9[) . using again the Laplace method. 
On the real axis, iph is treated separately in its oscillatory region [— 1 + 5, 1 — S], 
5 > 0, and in the neighbourhood of the turning points ±1. Hence, the sta- 
tionary phase method leads to an asymptotic expansion in powers of h of the 
denominator in (|2.9p . Finally, the statement follows from the Bohr-Sommerfeld 
quantization rule for h n (see [7], Exercise 12.3) or Wcyl formula [5]. □ 

In the harmonic case k — 1 (Davies operator), the first term in (|2.5D yields 
the Davies-Kuijlaars theorem [3] : 

1 f I \ ( e l0 ' A \ 

lim -log||n n || = Cl (0)=4<p 1 , == =2Re/ 



n^+oon °" r \ v /2cos(6»/2) I \ ^2 cos(6>/2) 



where f(z) = log(z + V^ 2 — 1) — zv? 



3 Eigenfunctions and semigroups 

The following theorem has been proved in [2] in the case of complex Airy 
operator 9), and in [4] in the harmonic case (k — 1), as well as for A(2k, 9), 
k > 2, \9\ < 2. The proof actually extends to any operator A(2k,9) with 



Theorem 3.1. For any m = 1, 2k, k > 1, and any 9 satisfying ji.^[ ) ; the 
eigenfunctions of A(m,9) form a complete set of the space L 2 (M). 

Notice however that the eigenfunctions of A(l,9) and A(2k,9), k > 1, can 
not form a Riesz basis because of the growth of the instability indices as n — > 
+oo. 

Theorem 13.11 and the previous estimates enable us to study the convergence of 
the operator series defining the semigroup e ~ t ^ m fi) associated with A(m,9) 
when decomposed along the projections II„. 

The following statement extends the result of [5] in the harmonic case. 
Corollary 3.2. Let \9\ < ir/2 and e -tA(m,e) be the 

semigroup generated by 

A(m,9), \ n ~\ n {m,9) the eigenvalues of A{m,9), and Yl n = Tl n {m,9) the 
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associated spectral projections. 

Let T{9) = ci(9) / cos(0/2), where c\{9) is the constant in h2.b)) . The series 



n=l 

is not normally convergent in cases m = 1 /or any t > 0, and m = 2 /or 
t < T(0) ; in cases m = 2 /or t > T(6), and m = 2k for any t > 0, k > 2, i/ie 
series converges normally towards e - t - 4 ( m . 6 ') arl ^ y or TV sufficiently large and 
for some constants C\ = C\{k, 9) and Ci = C2{9), the following estimate holds : 

11 1 <JVj " " 1 ^cxp(-2cos(0/2)(t-T(0))7V), k = l,t>T 

(3.1) 

where H < n = Hi + • • • + TIn-i denote the projection on the first N — 1 eigens- 
paces. 
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